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Abstract. We investigate the optical conductivity of a Hubbard ring in presence 
of an impurity by means of exact diagonalization using the Lanczos algorithm. We 
concentrate thereby on the first excited, open shell state, i.e. on twisted boundary 
conditions. In the metallic phase a substantial part of the spectral weight lies in 
the Drude peak, <y{u>) = D c S(lo) + cr reg ■ In the non-interacting system, the Drudc 
peak can be visualized in our calculations at lu = even for finite chain lengths. 
Adding the impurity, the main peak is shifted to finite frequencies proportional to 
the impurity strength. The shift indicates the energy gap of the disturbed finite 
size system, also in the interacting system. Thus, we can pursue in the optical 
conductivity for finite metallic systems the energy gap. However, due to level 
crossing, the impurity-induced peak arises in the interacting system first when a 
certain impurity strength is exceeded. In the Mott insulating phase, the impurity 
leads to an impurity state within the gap. 



1. Introduction 

The optical properties of strongly interacting electron systems are of great interest, 
primarily to understand the properties of metal-insulator transitions We 
consider here the optical properties of Hubbard rings. The attention is directed 
to onc-dimcnsional systems, because the investigation of low dimensional systems 
has already provided important insights E] • One-dimensional systems are studied 
extensively both experimentally and theoretically in the context of correlated systems 

In addition, they are of great interest in the context of disordered systems UJ. 
Furthermore, powerful theoretical tools exist for the theoretical treatment of one- 
dimensional systems, like bosonization |7j, Bethe-ansatz [Hj, conformal field theory 
0, and the density matrix renormalization group [U3|- Also exact diagonalization is 
more powerful in low dimensions since larger systems can be investigated. 

In order to study the localization due to disorder in the one-dimensional 
Hubbard-model with regard to a metal-insulator transition, we investigate the optical 
conductivity of a Hubbard ring with impurity. In our work we ask the following 
questions: How do the different energy scales - vp/L, e, A - interplay in the metallic 
phases? Here vf is the Fermi velocity, L the system size, e the characteristic local 
impurity strength, and A the energy gap. How does the energy gap in the insulating 
phase, and the absorption edge seen in the optical conductivity change with the 
impurity? Especially, we concentrate on a weak impurity and weak interaction. The 
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transition from a Mott-insulator to an Anderson-insulator was studied previously in 
case of strong disorder and strong interaction ^2 ^] . 

The outline of the article is as follows. First we discuss the known properties 
of the disordered Hubbard-model. Then we describe our numerical approach to the 
optical conductivity. In section 4 we present our numerical results. 



2. The Hubbard model 

2.1. The clean model 

The Hubbard model is the generally accepted prototypical model to describe 
the interplay between kinetic energy (— ► derealization) and local interaction (—> 
localization) for electronic systems. The Hubbard chain is, in addition, exactly solvable 
by means of the Bethe ansatz ^3]- The Hamiltonian of the Hubbard model is given 

by 

N N 

flHubb = - t ^2( c i,<7 C i+i,*+ h - c) +Uj2 n iA n hl- (!) 

i,(7 i 

The chain length is given by L = Na, a is the lattice spacing, and n — N e /N is 
the particle density. In addition, we define rij = + n^i, rrii = — Uji, and the 
magnetization M — Ylii^t) = -^T ^ -^l- I n the clean case, the model shows three 
phases. Phase one, the Luttinger liquid phase, arises for U > and away from half 
filling. Spin- and charge-excitations are those of a Luttinger liquid. The point {7 = 0, 
n = 1 (non-interacting electrons and half filling) also belongs to this phase. Phase two, 
which is called spin gap phase, occurs for U < 0, where the spin-excitation spectrum 
has a gap and the low-lying charge-excitations can be described by those of a Luttinger 
liquid 7 . The last phase, the Mott insulating phase, occurs for U > and half 
filling, where the charge excitations have a gap and the spin-excitations are those of a 
Luttinger liquid. In the clean case, the Hubbard model shows spin-charge separation 
which is characteristic for Luttinger liquids ^JE3- A longer ranged interaction, even 
if it is the more generic case, leads to additional ground state phases 



2.2. Impurities in interacting Fermi systems 

We introduce disorder in our model by adding local potentials ej to the above 
Hamiltonian, Eq. Q), 

N 

H = Hnubb + ^ e '" 4 ■ ( 2 ) 

i 

The most studied example of impurity effects is an interacting system of spinless 
fermions in the presence of a local potential scatterer of strength e, i.e. e, = e<5i,o- 
The behavior of this system is well known. In this case, the free motion of the 
fermions inside the ring is restricted mainly due to the backscattering at the impurity, 
±fcp —¥ ^fkp. As discussed by Kane and Fisher ^7] the impurity strength scales to 
zero (becomes transparent) for an attractive interaction and scales to infinity (becomes 
completely reflective) for a repulsive interaction, according to the renormalization 
group equation 

de 

- = {1-K)e. (3) 
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The Luttinger parameter K depends on the interaction such that oo > K > 1 for 
attractive interaction, K = 1 in the non-interaction system, and 1/2 < K < 1 for 
repulsive interaction. In the charge ordered insulating phase at half filling, which 
corresponds to the Mott phase of the Hubbard model, K = 0. 

In a Hubbard chain the impurity is accordingly known to be relevant in the 
Luttinger liquid phase. The scaling of the impurity |17j . 

^ = (2-^-^)6, (4) 

contains now the Luttinger parameters K c and K s for the charge and spin degrees of 
freedom. As before, K c = 1 and K s = 1 in the non-interacting system, K c > 1 for 
attractive interaction and K c < 1 for repulsive interaction. K s = in the spin gap 
phase K s — and K s = 1 for repulsive interaction. In contrast to spinless fermions, 
no metal-insulator transition as function of the interaction is found in the Hubbard- 
model in presence of the impurity; the system is always localized. The scaling beyond 
Eq. (J3J) and @ is discussed by Meden et al. JH| They show that the scaling of the 
impurity - as shown in Eq. Q and predicted by bosonization, is only valid for long 
chains and strong interaction. 

These results for a single impurity can be extended on disordered systems 6 . 
Thereby a metal-insulator transition is found in the spinless Fermi model at finite 
disorder strength for strong attractive interaction In the Hubbard model it is 

found that a weak repulsive interaction reduces the localizing effect of the disorder j5U] . 
However, it is also known [2] that the 4fcp contribution of the scattering potential, not 
included in the above scaling relation, is responsible for the localization in interacting 
systems. The Akp contribution becomes relevant for strong interaction. Whereas the 
2Uf impurity-scattering and interaction weaken each other, thus leading to weaker 
localization than in the non- interacting system, the 4fc^ contribution enhances the 
localization for strong interaction. 

Many recent activities related to the disordered Hubbard model concentrate on 
single defects. For example, the Friedel oscillations are investigated |ffl 1221 12*3*] . The 
exponent of the algebraic decay depends strongly on interaction and resembles the 
scaling of the impurity for a given interaction. Thus, studying the local behavior, 
possibly the relevant theoretical model can be identified [21] ■ Concerning the insulator- 
insulator transitions which occur in the disordered interacting system (Anderson 
insulator versus Mott insulator) [111 12*5] and the Peierls-Hubbard model (band 
insulator versus Mott insulator) O EH] , the Drudc weight - or the phase sensitivity - 
and the optical conductivity were used to characterize the different insulating phases. 

3. Optical conductivity 

In the linear response approximation, the optical conductivity gives the response of a 
system to an external electric field, (j) = crE, E = — A/c. The Kubo formula is used 
to evaluate the optical conductivity a(t, t') by means of the current-current correlation 
function < [j(f),j(*')] >, 

<j(*)> + Jdia = i / df <[j(t),j(f )]> - — A(t) , (5) 

J_ 00 c mc 
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with j = ep/m. Using the Green's functions formalism to evaluate the propagation of 
the electrons, the optical conductivity is given by 

— -i- ' ' 



Recr(w) 





E 



oj + id-H iq 



,\n)\ 2 8(Lu-(E n -E )) 



The current operator for a one-dimensional lattice model is given by 



(6) 
(7) 

(8) 
(9) 



j q = y^ jiexp(iql) = 2iexp(ig/2)^sin(fc + g/2)c£ CT Cfc_ gi0 . . 

I ka 

Using the fact that the Fermi surface of a chain consists of only the two Fermi points 
at ±fcp, the current operator can be written in an intuitive way, assuming linear 
dispersion relation near the Fermi points in the non-interacting system, e(k) oc vpk, 
vp = —It sm(kp ), 



j = v F 22(n R ^ k - n Ljk ) 



(10) 



where L, R denote left and right moving particles. Since j commutes with the kinetic 
and interaction part of the Hubbard Hamiltonian, it is a good quantum number. Using 
Eq. (JHJ the optical conductivity, Eqs. JSJ and (Q) can be rewritten as 



Imer (uj) = — 



kin 



2(0|j> 



{E n — Eq) 



(11) 



Recr(w) = DS(lu) + cr rcg . (12) 

The Drude weight D is given by D ~ 2v c K c in the Hubbard model, where v c and K c 
are interaction dependent constants which can be determined by means of the Bethe 
ansatz [Z7\ . The Drude weight is also related to the charge stiffness and the phase 
sensitivity |58j. Using Eq. (fTU|) . first the conjecture of persistent currents and infinite 
optical conductivity was given [23]. Finally, the optical conductivity obeys a sum rule, 

r°° ire 2 

/ Rea(w)dw = -(H kin ) . (13) 

Joo L 

Within the Lanczos algorithm, also called Lanczos vector method ^31 , the 
coefficients arising in the representation of the correlation function as a continued 
fractions can be determined in an exact diagonalization calculation. Extensions of this 
method to longer chains using the density matrix renormalization group treatment are 
discussed in [321 an d |3S1- They also give a comprehensive overview of the Lanczos 
algorithm. The Green's function in Eq. © can be rewritten in continued fractions: 

(0\j+j q \0) 



a(uj + id) 



uj + iS — oq — 



(14) 



id — a\ — 



The coefficients a n and b n can be evaluated from fi n = {j+ H n j q ). Using a projective 
technique |34| the coefficients are determined by the following recursion formulas: 

\fo) =i|0) (15) 
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|/„ +1 > = (H - a n )\f n ) - bUf,^} (16) 

a n = (fn\j\fn)/(fn\fn) (17) 

K = (fn\j\fn)/{fn-l\fn-l) ■ (18) 



4. Numerical results 

In the following, we present our numerical results, obtained for a half- or third-filled 
band, in the four different sectors of the Hubbard model. Furthermore, we assume 
M = 0, i.e. N± — Np We compare our results from the optical conductivity with the 
energy gap, 

A(L) = [E(L, N e + 1) + E(L, N e - 1) - 2E(N e )}/2. (19) 

In the non-interacting case we can visualize the Drude peak at oj — even for finite 
system size. In case of N e /2 odd and anti-periodic boundary conditions (cl+i = — c\) 
or N e /2 even and periodic boundary conditions {cl+i — c\) four degenerate states 
contribute to the ground state. Due to this degeneracy of the partly filled and open 
shell states dissipationless excitations exist. Nevertheless, the true - i.e. the state with 
the lowest energy - ground state of a finite size system is obtained for N e /2 odd and 
periodic boundary conditions or vice versa. The different configurations are shown in 
Fig. ^ The upper left corner shows the partly filled shell state with a single particle 
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Figure 1. Single particle ground state energy, versus momentum, fc, for 

periodic - k = 2niz/L - and anti-periodic - k = (2n + l)ir/L - boundary 
conditions. We compare the ground state with ten and twelve electrons on ten 
sites. The true ground state for ten electrons is plotted in the upper right corner, 
the true ground state for twelve electrons in the lower right corner. The spins 
are only drawn for the levels near the Fermi energy, and the plus and the small 
square correspond to f J.. For the partly filled shell state - upper left corner - only 
one of the spin-degenerate states is shown. For the open shell states - lower left 
corner - only one realization is shown. As indicated in the upper right corner, 
the line gives the dispersion, e = — 2tcosk of the free tight binding model in the 
thermodynamic limit. 
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at the Fermi points, the lower left corner shows the open shell state with two particles 
at one of the Fermi points, the upper right corner shows the filled shell state according 
to the true ground state for N e /2 odd, and the lower right corner the filled shell state 
according to the true ground state for N e /2 even. The statement that the closed shell 
state is always the lowest in energy holds also in the interacting model |2H], an d in 
the disturbed model |29| . 

The twist in the boundary conditions is often modeled by a magnetic flux inside 
the ring. Anti-periodic boundary conditions corresponds to cl+i = exp(i7r)ci 
Using this formulation, it is clear that the magnetic flux can induce a current in 
the system. The open shell states, which are realized for N e /2 odd and anti-periodic 
boundary conditions, see the plot in the bottom left corner of Fig. carry momentum 
and current, see Eq. (|10J) and jSH]. Therefore, they have an overlap with the current 
operator. In the following we show only results for boundary conditions according to 
the open shell states. 



4.1. U = 

As seen in Fig. [21 the Drude peak is found at u> = in the clean system if we choose 
anti-periodic boundary conditions. According to D — Trvp, the peak is higher in the 
half filled case than in the third filled case. For the true ground state with the filled 
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Figure 2. Optical conductivity in arbitrary units versus frequency in units of 
the hopping rate t. System size and electron number are L = 10, N e = 10 
in the half-filled case and L = 9, N e = 6 in the third filled case. In the first 
case we compare different impurity strengths, in the latter case we compare 
the ground state (periodic boundary conditions) with the open shell state (anti- 
periodic boundary conditions). In the inset the higher excitations are shown for 
the third filled band. 



shell no Drude peak or another prominent peak at higher frequencies with uj « 1/L is 
seen. The first peak at u> « 1.2t corresponds to the energy of adding a particle above 
the filled shell. Adding the impurity the energy levels of the non-interacting system 
are partly shifted to higher energies. 

In case of four lattice sites, the energies can be evaluated analytically. The 
open shell state with e = splits off in e = and e w 2e/L for small e < t. 
This relation applies also to longer chains. The impurity does not change the 
structure of the electron distribution in coordinate space. The difference in the 
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reorganization of the electron distribution due to an impurity using periodic or 
anti-periodic boundary conditions can be seen - for longer chains numerically - in 
the local density. In the ground state, i.e. for the filled shell, the electrons are 
equally distributed. An additional impurity induces decaying Friedel oscillations. In 
contrast, the electron distribution shows oscillating behavior - already without the 
impurity - if the degenerate open shell state is considered. Adding the impurity, the 
electron density at the impurity site is reduced and we see the Friedel oscillations in 
addition. Looking at momentum space we see, that the lower level corresponds to 
the asymmetric superposition of the wave functions at the fc-points at the Fermi-level 
|0) = V*-7r/2 — V-V/2- The higher level, |1) corresponds to the symmetric superposition. 
The current operator connects both states, (0|j|l) ^ 0. Thus, adding the impurity 
the main peak shifts to higher frequencies, with wp oc 2e/ L. With increasing impurity 
strength, the peak-height decreases slightly, with nearly constant width, thus the 
peak loses weight. This weight occurs dominantly in the next peak around to « 1.2t. 
In the thermodynamic limit, L — > oc, the peak position shifts to zero frequency, 
up ~ 1/L -> 0. 

4-2. Luttinger phase 

In the following, we want to discuss the interacting system regarding a third filled 
band. In the numerical treatment we consider six electrons on nine lattice sites. 
With interaction the ground state cannot be viewed within in the band scheme. The 
degeneracy of the states, which have contributed to the ground state of the non- 
interacting system, is lifted. The splitting is proportional to U/L. The ground state 
with anti-periodic boundary conditions is - roughly speaking - given by the symmetric 
superposition of partly filled shell states where the spins are in spin-singlet states. The 
state of the next level has no overlap with the ground state via the current operator. 

However, adding the repulsive impurity to the interacting system, there is an 
interplay of interaction and impurity. Due to the interaction the charge distribution 
tends to be homogeneous, whereas the repulsive impurity repels the electron from one 
site. In the band picture, the partly filled shell states split off due to the impurity. If 
the level splitting is smaller than the splitting due to interaction, no rearrangement 
is expected. Increasing the impurity strength a state similar to the non-interacting 
system with impurity becomes more favorable. In this case, there is an overlap between 
the lowest states. 

In the numerical data, the Drude peak occurs abruptly by increasing the impurity 
strength for fixed interaction strength, see Fig. [3J The transition is very sharp. The 
characteristic impurity strength is proportional to the interaction, £ C (U) oc U/2, as 
indicated in the inset of Fig. [3] By comparing with the energy gap, Eq. i|19|) . 
we see that the Drude peak occurs when Al(U) becomes independent of impurity 
strength, or Ai(e) shows a dip at e c , see Fig. 0] Once raised, the peak position 
is proportional to the gap, ujp oc A([/, e). We expect the Drude peak occurring 
for 2e ~ U at up « 2e e g/L due to the discussion of the last subsection. With 
e cff = e(L/L ) 1 - {K ^ +Ks)/2 |H|, we find lo p w 2eL 1 ~ u /^ v ^ . We are thus able to 
follow in the optical conductivity of finite metallic systems the energy gap on the 
basis the Drude peak. 
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Figure 3. Optical conductivity in arbitrary units versus frequency of the 
Hubbard model at third filling. System size and electron number are L = 9 
and N e = 6. The inset shows the characteristic impurity strength e c , where the 
peak first arises, versus interaction U. 
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Figure 4. a) Energy gap versus interaction and fixed impurity strength, b) 
Energy gap versus impurity strength and fixed interaction. In both cases, L = 9 
and N e = 6, and we use anti-periodic boundary conditions. 
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4-3. Spin- gap phase 

Within the metallic phase for attractive interaction, where the spin excitations are 
frozen out, the optical conductivity shows a similar behavior as in the Luttingcr liquid 
phase. Now, the superposition of the open shell states is always the ground state. In 
contrast to the Luttinger liquid phase, the Drude peak is present for all interaction 
and impurity strengths (in case of a repulsive impurity). The numerical results are 
shown in Fig. EK). The peak position marks again the energy gap, wp ~ A(U,e,L). 
Both increase first with \U\, but decrease then with l/|f|, as shown in Fig. [SJd). 
With periodic boundary conditions the gap decreases with l/{/-behavior over the 
whole parameter regime. Due to the spin gap, the peak loses drastically weight with 
increasing interaction strength. 
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Figure 5. a) Optical conductivity in arbitrary units versus frequency. L = 10, 
N e = 10. b) Energy gap and peak position — taken from a) — versus interaction. 
The lines are connecting the data points. 



4-4 ■ Matt insulating phase 

In case of half filling, the interaction lifts the degeneracy of the open shell states and 
the partly filled shell states [37|- Again we refer to the four site system with four 
electrons to write down the analytical expressions. The fourfold degenerate open shell 
states with E = -At splits off in E = -4t + 3J7/4, E = -At + U, and E = -At + W/A, 
where E = J2 n £n - ^ ne secon d state with E = —At + U is doubly degenerate and is 
represented by the two spin-degenerate partly filled shell states shown in the upper 
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left corner of Fig. ^ The ground state with interaction is given by the asymmetric 
superposition of the open shell states. For this reason, we see in the numerical data 
for half filling a peak at utp w A(L) w U, compare We will call this peak 

'interaction' peak. In the insulator the question arises, whether the often supposed 
equality between the energy gap and the frequency of the interaction peak persists 
also in the system with impurity. 

With weak impurity, compare the data for U — 2t and t — tin Fig. El the position 
of this main peak shifts to lower frequencies as the energy gap becomes smaller with 
increasing impurity strength. An exception is the case of small interaction, where the 
gap decreases slightly with impurity strength, but the frequency of interaction peak 
increases. For stronger interaction, both decrease with impurity strength. Comparing 
the data taken from the optical conductivity with the energy gap, it generally applies 
that interaction frequency and energy gap agree only for weak impurities. In addition 
to the interaction peak another peaks grows at small frequencies due to the impurity. 
The larger the interaction strength, the smaller is the weight of the impurity induced 
low-frequency peak. On the other hand, the stronger the impurity, the smaller is the 
weight of the interaction peak. 
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Figure 6. Optical conductivity in arbitrary units versus frequency. We use a 
smoothing of S = 0.1. System size and electron number is L = 10 and N e = 10; 
we compare U = t and U = 2t and e = and e = t, respectively. 

In case of small interaction, compare U — t and e = t in Fig. El or equivalently 
for strong impurities e w U, the behavior is opposite. The spectrum of the optical 
conductivity is governed by the impurity. The most surprising result is that the 
frequency of the impurity peak decreases to zero, Wi mp — ► for e ~ U, see Fig. At 
the same time, the peak gains weight. Increasing the impurity strength further, only 
the impurity induced peak survives and the interaction peak disappears. Thus, within 
the short chains we can see that the impurity can shift the lower lying peak in the 
spectra to zero frequency, thereby indicating a metallic state. Nevertheless, the energy 
gap is finite in this case. From this point of view, the impurity is relevant in order 
to destroy the Mott insulating state. Increasing the impurity strength further, only 
the impurity induced peak survives and the interaction peak is absent. The impurity 
induced peak then shifts to higher frequencies by further increasing the impurity 
strength - comparable to the behavior in the Luttinger phase. 

The results in the Mott insulating phase clarify the results of our previous 
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Figure 7. Energy gap and peak positions, aJA.t>> f° r the interaction, and ojj mp , for 
the impurity peak, versus impurity strength. The data are taken from Fig. 151 for 
U = t. The data for the energy gap are obtained again for L = 10 and N e = 10 
and anti-periodic boundary conditions. The lines are connecting the data points. 

investigations, where we found a fast decay of the Friedel oscillations for small 
interaction |23) . A similar behavior is known to exist for an irrelevant impurity in 
a metal. The results concerning the phase sensitivity, which show a slight maximum 
for small interaction in presence of a small impurity, indicates also a tendency to 
derealization. Due to the strong influence of the impurity in the small gap regime, 
the correlation gap is closed and gapless excitations are possible. The resulting 'metal' 
has other properties than the Luttinger liquid and we can find fast decaying Friedel 
oscillations. 



5. Summary and conclusions 

In summary, a detailed determination of the ground state characteristics of the one- 
dimensional Hubbard model in the presence of an impurity has been achieved. In 
particular, we have analyzed the optical conductivity of a small ring in all phases of 
the one-dimensional Hubbard model. In the non-interacting system the Drude peak is 
seen in small rings if appropriate boundary conditions are chosen, for which the ground 
state is degenerate. Nevertheless, the continuation to the thermodynamic limit, where 
the transport properties should be independent from the boundary conditions, is not 
clear. Adding the impurity to the non-interacting system, the degeneracy of the 
ground state is lifted and the optical conductivity shows a clear peak at ujp ~ e. 
Thus, so to speak, the peak position maps the behavior of the finite size energy gap. 
A similar behavior is found in the Luttinger phase and the spin gap phase of the 
interacting system. 

An impurity in the Mott insulating phase leads to an impurity state within the 
correlation gap. The qualitative behavior of the optical conductivity is different 
for small and large - compared to the interaction - impurity strength. For large 
interaction the impurity has only a weak influence on the system, the interaction peak 
is only slightly shifted and has clearly more weight than the impurity peak. On the 
other hand, in case of a small energy gap, hence for small interaction, the impurity 
dominates the behavior of a(to). The peak related to the impurity is only seen provided 
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e > U. As in previous studies, we find that interaction and impurity - both in the 
Luttinger and in the Mott phase - weaken each other. 
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